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Abstract 



Comparisons of cosmological models to current data show that the presence of a non- 
CS| , trivial feature in the primordial power spectrum of fluctuations, around the scale 

_ k ~ 0.05/iMpc~^, is an open and exciting possibility, testable in a near future. This 

^ ■ could set new constraints on inflationary models. In particular, current data favour 

. a ACDM model with a steplike spectrum, and more power on small scales. So far, 

this possibility has been implemented only in toy models of inflation. In this work, 
we propose a supersymmetric model with two U{1) gauge symmetries, associated 
with two Fayet-Iliopoulos terms. Partial cancellation of the Fayet-Iliopoulos by one 
of the charged fields generates a step in the primordial power spectrum of adiabatic 
perturbations. We show that when this field if charged under both symmetries, the 



OO 

^ • spectrum may have more power on small scales. 



1 Introduction 



A most exciting aspect of inflationary models is the possibility of constrain- 
ing them through observations of cosmological perturbations. Generally, the 
power spectrum of primordial fluctuations can be easily related to the infla- 
tion potential [1]. In more complicated models, like multiple-stage inflation 
[2-4], slow-roll conditions may be violated; then, the primordial spectrum has 
to be derived from first principles, and exhibits a non-trivial feature at some 
scale. These scenarios shouldn't be regarded as unrealistic: inflation has to 
take place within the framework of a high-energy theory, with many scalar 
fields. So, the real question is not to know wether multiple- stage inflation can 
occur, but wether it is predictive, or just speculative. Multiple-stage inflation 
has already been invoked for many purposes : generation of features in the 
power spectrum at observable scales, setting of initial conditions, description 
of the end of inflation. Here we are interested in the first possibility. The 
range of cosmological perturbations observable today in the linear or quasi- 
linear regime has been generated during approximately 8 e-folds, 50 e-folds 
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before the end of inflation. So, a feature in the primordial power spectrum 
could very well be observable, without any fine-tuning; only observations can 
rule out this possibility. 

In the radiation and matter dominated Universe, primordial fluctuations of the 
metric perturbations couple with all components, leading to the formation of 
Cosmic Microwave Background (CMB) anisotropics and Large Scale Structure 
(LSS). It is well established that a pure Cold Dark Matter (CDM) scenario, 
in a flat universe, and with a scale-invariant or even a tilted primordial spec- 
trum, is in conflict with LSS data. So, many variants of this scenario have 
been considered, including a Broken Scale Invariant (BSI) power spectrum. 
Indeed, standard CDM predictions arc improved when less power is predicted 
on small scales. Specific cases have been compared accurately with both LSS 
data and recent CMB experiments, including the double- inflationary spectrum 
of [5], and the steplike spectrum of [6]. However, even with such spectra, it has 
been shown [7-9] that standard CDM cannot be made compatible with all ob- 
servations. Independently, recent constraints from la supernovae [10] strongly 
favour the ACDM scenario, with a cosmological constant A ~ 0.6 — 0.7 and 
a Hubble parameter h = 0.6 — 0.7. In this framework, a reasonable fit to all 
CMB and LSS data can be obtained with a flat or a slighly tilted spectrum. 
So, ACDM is very promising and should be probably considered as a current 
standard picture, from which possible deviations can be tested. In this respect, 
in spite of large error bars, some data indicate a possible sharp feature in the 
primordial power spectrum around the scale k ~ 0.05/iMpc^^. First, APM 
redshift survey data seem to point out the existence of a step in the present 
matter power spectrum [11,3]. Second, Abell galaxy cluster data exhibit a sim- 
ilar feature [12] (at first sigh, it looks more like a spike than like a step, but in 
fact a stcplikc primordial spectrum multiplied with the ACDM transfer func- 
tion reproduces this shape [8]). Third, this scale corresponds to / ~ 300, i.e. to 
the first accoustic peak in the CMB anisotropies, and increasing power at this 
scale, through a bump or a step, would lead to a better agreement with Saska- 
toon. In the next years, future observations will either rule out this possibility 
(which will then be attributed to underestimated errorbars), or confirm the 
existence of a feature, and precisely constraint its shape. In [8,9], we compared 
a specific BSI ACDM model with CMB and LSS data. The primordial spec- 
trum was taken from a toy model proposed by Starobinsky [6], in which the 
derivative of the inflaton potential changes very rapidly in a narrow region. In 
this case, the primordial spectrum consists in two approximately flat plateaus, 
separated by a step. This model improves the flts to CMB and LSS data (even 
when the cluster data are not taken into account). It also reproduces fairly 
well the feature of Einasto et al. [12] when the step is inverted with respect 
to the usual picture, i.e., with more power on small scales. Independently, in 
a preliminary work [13] (motivated by the inflationary model of [14]), a pri- 
mordial spectrum with a bump centered at /c ~ 0.06/iMpc~^ was compared 
with CMB and LSS data, including more redshift surveys than in [8] , and not 
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taking Einasto et al. [12] data into account. It is remarkable that among many 
cosmological scenarios, this BSI spectrum combined with ACDM yields one 
of the best fits. 

So, we have good motivations for searching inflationary models based on realis- 
tic high-energy physics, that predict a bump or an inverted step at some scale, 
and approximately scale- invariant plateaus far from it. Successfull comparison 
with the data requires the deviation from scale invariance to be concentrated 
in a narrow region ki < k < k2; roughly, /c2/^i should be in the range 2-10. 
This is quite challenging, because double inflationary models studied so far 
predict systematically less power on small scale, with a logarithmic decrease 
on large scale rather than a plateau. However, these models were based on the 
general framework of chaotic inflation. Today, the best theoretically motivated 
framework is hybrid inflation [15]. Indeed, hybrid inflation has many attrac- 
tive properties, and appears naturally in supersymmetric models: the inflaton 
field(s) follow(s) one of the flat directions of the potential, and the approx- 
imately constant potential energy density is provided by the susy-breaking 
F or D-term. When the F-term does not vanish, conditions for inflation are 
generally spoiled by supergravity corrections, unless interesting but rather 
complicated variants are considered (for a very good review, see [16]). On 
the contrary, the simple D-term inflation mechanism proposed by Binetruy & 
Dvali [17] and Haylo [18] is naturally protected against supergravity correc- 
tions, and can be easily implemented in realistic particle physics models, like 
the GUT model of reference [19], without invoking ad hoc additional sectors. 
If supergravity is considered as an effective theory from superstrings, D-term 
inflation is hardly compatible with heterotic strings theory [16,20], but con- 
sistent with type I string theory [21]. 

Our goal is to show that in this framework, a simple mechanism can generate 
a steplike power spectrum with more power on small scales. This mechanism 
is based on a variant of D-term inflation, with two Fayet-Iliopoulos terms. 
However, the fact that it is D-term inflation is not crucial for our purpose: a 
similar lagrangian could be obtained with two non-vanishing F-terms, or one 
D-term plus one F-term, like in [4] . We will not consider the link between the 
model of this paper and string theory, putting by hand the Fayet-Iliopoulos 
terms from the beginning. In a very interesting paper, Tetradis & Sakellariadou 
[4] studied a supersymmetric double inflationary model with a quite similar 
lagrangian. However, the motivation of these authors is to save standard CDM. 
So, they are pushed to regions in parameter space quite different from us, and 
do not consider a steplike spectrum with flat plateaus, but rather a power-law 
spectrum with n = 1 on large scales and n < 1 on small scales. 
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2 The model 



We consider a supersymmetric model with two gauge symmetries U{1)a x 
U{1)b- and two associated Fayet-lliopoulos positive terms a-nd (there is 
no motivation from string theory to do so, at least at the moment, but SUSY 
and SUGRA allow an arbitrary number of Fayet-lliopoulos terms to be put by 
hand in the lagrangian). In this framework, the most simple workable model 
involves six scalar fields: two singlets A and B, and four charged fields A±, B±, 
with charges (±1,0), (±1,±1). Let us comment this particular choice. First, 
the presence of two singlets is crucial. With only one singlet coupling to both 
A-t and 5-J-, we would still have double-inflation, but the second stage would 
be driven by both F and D-terms, and no sharp feature would be predicted in 
the primordial spectrum. Second, each charged field could be charged under 
one symmetry only; then, a steplike spectrum would be generated, but with 
necessarily less power on small scales. Here, the fact that B_ has a charge 
— 1 under both symmetries is directly responsible for the inverted step, as 
will become clear later. Finally, both global susy and supergravity versions of 
this model can be studied: supergravity corrections would change the details 
of the scenario described thereafter, but not its main features. We consider 
the superpotcntial: W = aAA^A_ + l3BB^B_ (with a redefinition of A and 
B, we have suppressed terms in BA^A_, AB^B_, and made {a, (3) real and 
positive). In global susy, the corresponding tree-level scalar potential is: 



y=Q;2|^|2(|A+|?f|^_|2) + Q;2|A+^_|2 + M(|A+|?-|A_|2+|S+|2-|S_|VeA)' 

+/32|5|2(l5+lV|5-r)+/32|5+5-r+^(l5+r-li^-r+eB)^ (1) 

with a global supersymmetric vacuum in which all fields are at the origin, 
excepted = -\/Cb, and, depending on the sign of (^a — Cb); I^-I oi" I^+I — 



3 The two slow-roll inflationary stages 



3.1 Inflationary effective potential 



There will be generically two stages of infiation, provided that initial condi- 
tions for [A, B) obey to : 
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Then, charged fields have positive squared masses and stand in their (local) 
minimum A± = B± = (for a discussion of the charged fields initial condi- 
tions, see for instance [22,23]). A and B have a constant phase, while their 
moduli A = \A\/\/2 and B = \B\/\/2 behave hke two real infiaton fields and 
roll to the origin, until one inequality in (2) breaks down. We assume that the 
condition on B breaks first. 

During this first stage, the field evolution is driven by the effective potential: 
Vi = {g\il+g%^l)/2+AVi. The one-loop correction AVi is small (A^i < Vi) 
[21], but crucial for the field evolution. It consists in two terms with a logarith- 
mic dependence on A and B. The former takes a simple form following from 
[24], because we can assume g\C,A ^ The latter is more complicated 

because the dimensionless quantity b = {QaCa + 9b^b) / {P^lB]"^) goes to one 
when B reaches its critical value. The exact expressions are: 



^^^=i^(^"4^+2j 

327r2 V 2fe2 ) 

(A is the renormalization energy scale at which qa and qb must be evaluated). 
When 6^1, the term involving h tends to |. Even at 6 = 1, this term only 
increases the derivative [dVi/dB) by a factor 21n2 : so, it can be neglected in 
qualitative studies. In this approximation, it is easy to calculate the trajecto- 
ries of A and S, and to note that B reaches its critical value before A only if 
the initial field values obey to: 



This condition is natural, in the sense that it allows \A\q and |i?|o to be of the 
same order of magnitude, whatever the values of gA,B and $,a,b- 

At the end of the first stage, (5, -B_) evolve to another false vacuum: i? = 0, 
= {g\iA + 9B^B)/{gA + 9b)- During this transition, the charged fields 
B+, A± remain automatically stable if we impose ^b < 2^^. Afterwards, a 
second stage occurs: A rolls to the origin, driven by the potential: 




until or \A_\ becomes unstable, and quickly drives the fields to the su- 
persymmetric minimum. 
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3.2 Second stage of single-field inflation 



Let us focus first on tfie second stage of inflation, in order to find the small 
scale primordial power spectrum {i.e., if the second stage starts at t = t2, and 
k2 = a(t2)H{t2), the power spectrum at scales k > k2). This stage should last 
approximately ~ 50 e-folds, so that the transition takes place when scales 
observable today cross the Hubble radius. 

A standard calculation shows that for a of order one, the second slow-roll con- 
dition breaks before A reaches its critical value (which is given by a^l^cp = 

-t^-t"I'Ca ~ ^b\)- Integrating back in time, we find that N e-folds before the 
end of inflation, 

1^1 = J^S£Es^Up (6) 



(we are using the reduced Planck mass Mp = (SttG)"^/^ = 2.4 x lO^^GeV). 
Then, the primordial spectrum can be easily calculated, using the single-field 
slow-roll expression [1]: 

" TSTT^Mf, {dV2ldAy ~ 75M% g\gl ^ ~ ' ^ ^ 



To normalize precisely this spectrum (7) to COBE, it would be necessary to 
calculate the contribution of cosmic strings generated by symmetry breaking 
[25] , to make assumptions about the geometry and matter content of the uni- 
verse, and to fix the amplitude of the step in the primordial power spectrum 
(between COBE scale and ^2)- However, if perturbations generated by cosmic 
strings are not dominant, and if the primordial spectrum is approximately 
scale-invariant as required by observation, we can estimate the order of magni- 
tude of the primordial power spectrum at all observable scales: ~ 5 x 10~^°. 
So, at A; = A;2, inserting (6) into (7), we find the constraint: 

sJ\iA-iB\ - 3 X 10-^Mp ~ 10^^-^^GeV. (8) 



At first sight, this constraint could be satisfied when (\/^, \/Cb) ^-^e both 
much greater than 10~^Mp, and very close to each other. Then, however, the 
amplitude of the large-scale plateau would violate the COBE normalization, 
as can be seen from the following. Also, there is no reason to believe that one 
term is much smaller than the other: this would raise a fine-tuning problem. 
So, we will go on assuming that both Fayet-lliopoulos terms have an order of 
magnitude ~ \/Cb ~ 10~^Mp, just as in single D-term inflation. 
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The spectrum tilt 715 at A; = /c2 can be deduced from the slow-roll parameters 
(e, 7]) [1]. Like in any model of single-field D-term inflation, e can be neglected, 
and ng{k2) = 1 + 277(^2) = 1 — '^/N ~ 0.98. So, the spectrum on small scales 
is approximately scale-invariant. 



3.3 First stage of two-field inflation 



During the first inflation, the primordial spectrum calculation must be done 
carefully. If slow-roll conditions were to hold during the transition between 
both inflationary stages, the evolution of metric perturbations (for modes out- 
side the Hubble radius) would be described at first order by the well-known 
slow-roll solution (see for instance [26]): 

( dAVA + dsVB \ . . 

'''h^ ""dty Va + Vb J' ^' 



where $ is the gravitational potential in the longitudinal gauge. Here, Ci is 
the time-independent coefficient of the growing adiabatic mode, while dA and 
ds are coefficients related to the non-decaying isocurvature mode (in fact, 
only dA — ds is meaningful [26]). The expression of Va and Vb at a given 
time can be calculated only if the whole field evolution is known, between the 
first Hubble crossing and the end of infiation. Formally, in the general case 
of multiple fields 0j, the V^'s can be found by integrating dVi = {dV / d4>i)d4>i 
back in time, starting from the end of inflation. This just means here that 
during the second slow-roll, Va — V2, Vb — 0, and during the flrst slow-roll: 




'^{91 + 9%) 

2(,i + ,|) + 32^^ r^ + 2j- ^''^ 

We see immediately from (9) with Vb = that the isocurvature mode is 
suppressed during the second inflationary stage. On the other hand, it must 
be taken into account during the flrst stage. This leads to the well-known 
expression [26]: 

5^ (k) = ^ ( VI \ 

^ lhT^''^p\{dVAldAf {dVBldBfj^: ^ ^ 



where the subscript k indicates that quantities are evaluated at the time of 
Hubble radius crossing during the flrst stage. 
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If slow-roll is briefly disriipted during the transition (this is the interesting 
case if we want to generate a narrow step or bump), the solution (9) doesn't 
hold at any time, but we still have a more general exact solution, describing 
the adiabatic mode (in the long-wavelenght regime): 



Ci^-^ j adt^ + (other modes). (12) 



If, during the transition, the other modes are not dominant (as usually ex- 
pected), the matching of the three solutions : (9) before the transition, (12) 
during the transtion, and (9) afterwards, shows that Ci is really the same 
number during all stages: the slow-roll disruption doesn't leave a signature 
on the large-scale power spectrum (11). On the other hand, if a specific phe- 
nomenon amplifies significantly isocurvature modes during the transition, the 
same matching shows that during the second stage, there will be an addi- 
tional term contributing to the adiabatic mode. This role may be played by 
parametric amplification of metric perturbations, caused by oscillations of 
B^. Generally speaking, the possibility that parametric resonance could affect 
modes well outside the Hubble radius is still unclear, and might be important 
in multi- field inflation [27]. In our case, this problem would require a care- 
full numerical integration, and would crucially depend on the details of the 
one-loop effective potential during the transition. So, we leave this issue for 
another publication, and go on under the standard assumption that expression 
(11) can be applied to any mode that is well outside the Hubble radius during 
the transition. 

Let us apply this assumption to our model, and flnd the large scale primordial 
power spectrum {i.e., if the flrst stage ends at t = ti, and ki = a(ti)H(ti), 
the power spectrum at scales k < ki). The contribution to 6h arising from 
perturbations in A reads: 

_ V VI 9%{9ieA+9ieB){u-iB)\ , 



The transition lasts approximately 1 e-fold (indeed, during this stage, the evo- 
lution is governed by second-order differential equations with damping terms 
+3H(j)i, and {B,B-) stabilize within a time-scale At ~ H~^). During that 
time, A is still in slow-roll and remains approximately constant. So, using eqs. 
(7) and (13), it is straightforward to estimate the amplitude of the step in the 
primordial spectrum, under the assumption that 5H{k2)\B is negligible: 

„2 _ suh)\A (1 - ^B/uni + 9ieB/9ieA) 
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Since we already imposed < '^^a, P can easily be smaller than one, so that 
we obtain, as desired, more power on small scales. The simple explanation 
is that with B_ charged under both symmetries, the transition affects not 
only the dynamics of {B,B_), but also the one-loop correction proportional 
to (ln|A|^), in such way that the slope (dV/dA) can decrease by the above 
factor p. 

However, perturbations in B must also be taken into account. Their contribu- 
tion to the large scale primordial spectrum reads: 

x2| - y Vb 16vr M^^ia + ^7ia) |p|2 .... 



75T:m% {dVB/dBy 75M% {g\ + gl) 



At the end of the first stage, the value of \B\ is roughly given by M^p{d'^Vi/dB'^) — 
Vi, since after the breaking of this slow-roll condition, eq. (11) is not valid any- 
more, and B quickly rolls to its critical value. This yields: 

\B\ = sUa + 9%Cb ^^g^ 



We can now compare Sfj\A and 5|^|_b at the end of the first stage. It appears 
that for a natural choice of free parameters {qa ~ gs, ~ ^b ^ Ha — 
Cb\), both terms can give a dominant contribution to the global primordial 
spectrum. If Sjj\B dominates, we expect p > 1, and a tilted large-scale power 
spectrum (due to the decrease of which is at the end of its slow-roll stage). 
On the other hand, if Sfjls is negligible, the step amplitude is directly given 
by (14), and the large scale plateau is approximately flat, with ns — 0.98 like 
in single-field D-term infiation. As we said in the introduction, this latter case 
is the most interesting one in the framework of ACDM models. A numerical 
study shows that S^Ia ^ ^h\b holds in a wide region in parameter space. 
Indeed, we explored systematically the region in which 0.1 < gs/gA 10, and 
iVU: \/Cb) are in the range (0.1-10) x lO'^Mp. We find that > 10 ^h\b 
whenever: 

(^B > O.S^A, a/cI > I-Ia/C^) or (ygA^gB, J^>l.lJu)- (17) 



So, inside these two regions, the primordial spectrum has two approximately 
scale- invariant plateaus {ns ^ 0.98), and the step amplitude is given by (14). 

Further, a good agreement with observations requires a small inverted step, 
0.75 < p < 0.85 [8], and of course, a correct order of magnitude for the 
amplitude, Sj^iki) ~ O(10^^°). These additional constraints single out two 
regions in parameter space: 
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(2.2gA <9B< ^9A, 13(v^ - 4.5 x 10"^) < 10^^^ < - 2.5 x lO"^)) 

or (^5b < 1.5C/A, 150^^-2.5x10"^ < 90^A- 4.3x10-3^ . (18) 



In this section, we only studied the primordial spectrum of adiabatic pertur- 
bations. Indeed, it is easy to show that tensor contributions to large-scale 
CMB anisotropics are negligible in this model, like in usual single-field D- 
tcrm infiation (a significant tensor contribution would require gA,B ~ 0(10) 
or greater, while the consitency of the underlying SUSY or SUGRA theory 
requires gA,B ~ O(10~^) or smaller [21]). 



3.4 Transition between slow-roll inflationary stages 



We will briefly discuss the issue of primordial spectrum calculation for ki < 
k < During this stage, A is still slow- rolling, but B and B_ obey to the 
second-order differential equations in global susy: 

+ 3if£(_) + -2^ = 0. (19) 



The derivatives of the potential are given by the tree-level expression (1), plus 
complicated one-loop corrections. In supergravity, the trcc-lcvcl potential will 
be different from (1), even with a minimal gauge kinetic function and Kahler 
potential. With a non-minimal Kahler potential, an additional factor will also 

multiply both inertial and damping terms -B(-) + 3HB(^_y The evolution of 
background quantities has been already studied in a similar situation in [22]. 

As far as perturbations are concerned, the simplest possibility would be to 
recover the generic primordial power spectrum introduced by Starobinsky [6], 
also in the case of a transition between two slow-roll inffationary stages, with a 
jump in the potential derivative. This would be possible if : (i) H was approx- 
imately constant in the vicinity of the transition; (ii) the number of e-folds 
separating both slow-roll regimes was much smaller than one (in other terms, 
the time-scale of the transition must be much smaller than H~^). However, 
in the model under consideration, (i) is not a good approximation, because 
during the transition there is a jump in the potential itself. Moreover, (ii) is 
in contradiction with the equations of motion (19) : B and S_ only stabilize 
after a time-scale 5t ~ ^ due to the damping factor +3iJ. This statement 
is very general, and holds even in supergravity with a non-minimal Kahler 
potential. So, it seems that only a first-order phase transition can reproduce 
the exact power spectrum of [6], as previously noticed by Starobinsky himself 
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[28] . A model of inflation with a first-order phase transition has been proposed 
in [14]. 

Further, we don't believe that the power spectrum in our model can be ap- 
proximated by the numerical solution of double chaotic inflation [5]. Indeed, 
in chaotic double inflation, fields masses are typically of the same order as the 
Hubble parameter: so, background fields have no strong oscillatory regime. 
In our case, estimating the effective mass for i?_ from the tree- level poten- 
tial, and comparing it with V (which gives a lower bound on H), we see that 
before stabilization, S_ will undergo approximately ~ Mp/-v/|s ~ 10^ oscilla- 
tions. These oscillations are relevant for the primordial spectrum calculation, 
because B_ perturbations and metric perturbations will strongly couple. 

So, for each particular model, a numerical integration of the background and 
perturbation equations should be performed in order to find the shape, and 
even the width of the transition in the primordial power spectrum. The result 
could be either a step or a bump, and since the fields stabilize in 5t ~ H~^, 
it is reasonable to think that the width of the feature will not be too large, 
(^2/^1) < 10. So, in any case, the primordial spectrum should be in good 
agreement with current observations, but precise comparison with future data 
requires numerical work. 



4 Conclusion 

We introduced a model of double D-tcrm infiation, with two f/(l) gauge sym- 
metries, and two associated Fayct-Iliopoulos terms. A phase of instability for 
two fields {B, B_) separates two slow-roll inflationary stages. During this tran- 
sition, S_ partially cancels the Fayet-Iliopoulos terms, causing a jump in the 
potential; also, since it is charged under both symmetries, it affects the one- 
loop corrections in such way that the potential can become less steep in the 
direction of one inflaton field, A. As a result, for a wide window in the space of 
free parameters (the Fayet-Iliopoulos terms ^a,b and the gauge coupling con- 
stants gA,B), the primordial spectrum of adiabatic perturbations consists in 
two approximately scale-invariant plateaus, separated by an unknown feature, 
presumably of small width, and with more power on small scales. The ampli- 
tude of the step, p, is given by a simple function of parameters (14). In the 
framework of ACDM, spectra with 0.75 < p < 0.85 are likely to fit very well 
current LSS and CMB data, as argued in [8,9]. However, for a detailed com- 
parison of our model with observations, we need the shape of the primordial 
power spectrum between the plateaus. This issue requires a numerical inte- 
gration, and the result will be model-dependent, in contrast with predictions 
for the slow-roll plateaus. Finally, before any precise comparison with observa- 
tions, one should also consider the production of local cosmic strings, which is 
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a typical feature of D-term inflation [25]. Indeed, CMB anisotropies and LSS 
may result from both local cosmic strings and inflationary perturbations [29]. 
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